Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



The red tail of carbon stars in the LMC: 
Models meet 2MASS and DENIS observations 

CO 1 P. Marigo 1 , L. Girardi 2 , & C. Chiosi 1 

O ' 

O - 

£SJ ■ Dipartimento di Astronomia, Universita di Padova, Vicolo dell'Osservatorio 2, 35122 Padova, Italy 

^ ' 2 Osservatorio Astronomico di Trieste, Via Tiepolo 11, 34131 Trieste, Italy 

Fjr^ ■ To appear in A&A, see http://pleiadi.pd.astro.it for a version with high-resolution figures. 

00 

Abstract. Carbon stars are known to exhibit systematically redder near-infrared colours with respect to M-type 
stars. In the near-infrared colour-magnitude diagrams provided by the 2MASS and DENIS surveys, the LMC 
C-type stars draw a striking "red tail", well separated from the sequences of O-rich giants. So far, this conspicuous 
feature has been absent from any set of available isochrones, even the few existing ones that include the TP-AGB 
evolution of low- and intermediate-mass stars. To investigate such issue we simulate the complete 2MASS K s vs. 
(J — K B ) data towards the LMC by means of a population synthesis approach, that relies on extended libraries 
of published stellar evolutionary tracks, including the TP-AGB phase. The simulations provide quite a detailed 
description of the several vertical "fingers" and inclined sequences seen in 2MASS data, due to both Galactic 
foreground and LMC O-rich stars. Instead, as mentioned, the red tail of C-stars sets a major difficulty: we find 
that TP-AGB models with solar-scaled molecular opacities, the usual assumption of existing AGB calculations, 
do not succeed in reproducing this feature. Our tests indicate that the main reason for this failure should not be 
• ascribed to empirical T e s-{J — K) transformations for C-type stars. Instead, the discrepancy is simply removed by 

Oh' adopting new evolutionary models that account for the changes in molecular opacities as AGB stars get enriched in 

carbon via the third dredge-up (Marigo 2002). In fact, simulations that adopt these models are able to reproduce, 
5_l ■ for the first time, the red tail of C-stars in near-infrared CMDs. Finally, we point out that these simulations 

also provide useful indications about the efficiency of the third dredge-up process, and the pulsation modes of 
long-period variables. 
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1. Introduction of a red plume of C-stars in a striking way (Cioni et al. 

1999; Nikolaev & Weinberg 2000). This feature can be 

It has been already known for a long time that carbon ap p re ciated in the K s vs. (J - K s ) diagrams of Fig. □ 

(C-type) stars occupy a particular region of near-infrared they show a marked a l m ost-vertical sequence of red gi- 

colour-colour diagrams (e.g. Richer et al. 1979; Cohen ants at (J _ K j _ Llj from which an inclined branch 

et al. 1981; Frogel et al. 1990; Hughes & Wood 1990). departs at Ks „ n towards redder colours, extending up 

Specifically, C-type stars exhibit systematically redder to (j_^ s ) _ 2 .0. This inclined branch is what we refer to 

colours than M-type stars, and particular colour-colour ag « the red tail „ Nikolaev & Weinberg (2000) and Cioni 

diagrams (e.g. J — H vs. H ~ K) display such a sharp et aL (20Q1) demonstrate that the red tail is populated 

dichotomy between the two spectral classes that they are by non _ obscured C stars, whereas the almost-vertical fea- 

often used diagnostic tools to identify carbon star candi- ture a t (J — K s ) ~ 1.1 is composed basically of O-rich 

dates from observed samples. giants ( of spectra i types from late-K to M). Moreover, the 

Recent wide-area near-infrared surveys like the Two luminosities of red tail stars indicate that they belong to 

Micron All Sky Survey (2MASS; Skrutskie et al. 1997) the thermally pulsing asymptotic giant branch (TP-AGB) 

and the Deep Near-Infrared Southern Sky Survey (DENIS; phase, and hence they are expected to be the result of car- 

Epchstein et al. 1997) have confirmed these findings. By bon surface enrichment by the recurrent third dredge-up 

providing photometric data for the complete AGB popu- process, 
lation of the Magellanic Clouds, they reveal the presence 

As for the modelling of this feature, the present-day 

Send offprint requests to: P. Marigo situation is as follows. Most of the stellar isochrones for 

e-mail: mar igo@pd . astro . it old and intermediate ages available in the literature do not 
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Fig. 1. Left panel: K vs. (J - K) diagram for the LMC region of sky (101° < a < 61° and -77° < S < -63°). The 
data are taken from the 2MASS Second Incremental Data Release. Note the red tail of visible C-stars departing from 
the sequence of red giants at (J — K) ~ 1.4 and extending up to (J — K) ~ 2. Right panel: The same from DENIS 
Point Source Catalogue towards the LMC (Cioni et al. 2000). 



include the TP-AGB phase, as they usually extend up to 
tip of the red giant branch (RGB) or the Early- AGB (e.g. 
Lejeune & Schaerer 2001; Dommguez et al. 1999; Weiss 
& Salaris 1999; Bergbusch & VandenBerg 2001 and refer- 
ences therein). Two main reasons for this discrepancy are 
that i) fully modelling the AGB phase is indeed complex 
and time-consuming, and ii) most stellar models fail to 
explain the formation of carbon stars except for those at 
higher luminosities. 

To our knowledge the only available isochrone sets that 
include TP-AGB stars are the Padova ones (Bertelli et a. 
1994; Girardi et al. 2000; Salasnich et al. 2001; Marigo & 
Girardi 2001). In these models the entire TP-AGB evo- 
lution is followed - with the aid of a synthetic approach 
to various degrees of detail - up to ejection of the stel- 
lar envelope. However, all of them fail to predict the red 
tail of C stars, and in fact the reddest TP-AGB tracks 
hardly reach (J — K) colours as large as 1.3. This is illus- 
trated in Fig. [5] If in the older isochrone sets (Bertelli et 
al. 1994; Girardi et al. 2000) the basic problem was the 
lack of the C-star phase in the stellar models, the same 
does not affect the most recent models (Marigo & Girardi 
2001), that account for the formation of C-stars via the 
third dredge-up. 

What is then the problem with the colours of theoret- 
ical C-star models? Why do they not reach (J — K) as 
red as ~ 2? Since in Marigo & Girardi (2001) isochrones 
the third dredge-up process is calibrated so that C-stars 
appear at the right luminosities (see Marigo et al. 1999 
for details), the problem might be related either to (i) an 
improper transformation from T e ff to the (J — K) colour, 



which introduces errors of several tenths of a magnitude 
for C-type stars, and/or (ii) an improper modelling of the 
C-star radii and effective temperatures, with typical errors 
of some hundreds of degrees Kelvin. 

In case of alternative (i) , a simple revision of the T e ff- 
colour relations, via e.g. the use of suitable relations for M- 
and C-type stars, would solve the problem. In the case of 
option (ii), the implications would be much deeper, since 
the T e ff values of Marigo et al. (1999) models are typical 
of most TP-AGB models in the literature. 

Exploring both possibilities is the main scope of this 
paper, which is organised as follows. In Sect.Glwe start by 
recalling the main prescriptions in our synthetic TP-AGB 
models, that also account for the transition from O-rich 
(M-type) to C-rich (C-type) stars (Sect. 12. l|) . A brief sum- 
mary of the present status of other AGB evolution models 
available in the literature is given in Sect. 12.21 putting par- 
ticular emphasis on the role of molecular opacities. 

Section[S]describes our simulations of the 2MASS data 
in the K vs. (J — K) colour magnitude diagram (CMD). If 
several 2MASS features are reproduced in detail, the red 
tail of C-stars is missed by our initial models. Section 0] 
investigates the possible reasons for it. Finally, our main 
conclusions are summarised in Sect. 

2. Models for C-stars 

2.1. Our prescriptions: old and new 

In this study the formation and evolution of C-stars is 
described with the aid the synthetic TP-AGB model de- 
veloped by Marigo et al. (1996, 1998, 1999), and Marigo 
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Fig. 2. A subset of theoretical isochrones described in 
Marigo & Girardi (2001), and transformed to JK pho- 
tometry by means of Girardi et al. (2002) tables. TP- 
AGB models are computed with the Kfi x prescription. The 
metallicity is Z = 0.008, and isochrones are presented at 
equally-spaced age intervals of Alogt = 0.1. The dotted 
lines mark the isochrone sections corresponding to C-stars 
(with C/0> 1). The plot limits correpond closely to those 
of Fig.Q]for a LMC distance modulus of 18.5 mag. In these 
and other models in literature, C stars do not form any 
red tail. 

(2002) to whom the reader is referred for all the details. 
For the sake of clarity let us briefly outline the main struc- 
ture of the model and its basic prescriptions. 

— The TP-AGB evolution - of a star of given initial mass 
and metallicity - is followed from the first thermal 
pulse up to the complete ejection of the envelope by 
stellar winds. This is performed by suitably coupling 

(i) analytic prescriptions (e.g. the core mass-luminosity 
relation, the core mass-interpulse relation, the growth 
rate of the core mass during the quiescent stages), and 

(ii) a complete static envelope model to be integrated 
from the atmosphere down to the core. 

— The initial conditions at the first thermal pulse 
are derived from sets of stellar tracks for low- and 
intermediate-mass stars by Girardi et al. (2000), which 
cover the previous evolutionary phases, i.e. from the 
zero- age main sequence up to the beginning of the TP- 
AGB phase. In this way consistency and homogeneity 
are guaranteed in the construction of the correspond- 
ing stellar isochrones (see Sect. I3.2)l . 

— Third dredge-up episodes possibly take place at ther- 
mal pulses and eventually lead to the transition from 
M- to C-class if the surface C/O ratio increases above 



1. The third dredge-up is described as a function of 
two parameters, namely: the efficiency A, and the min- 
imum temperature T^ e d at the base of the convec- 
tive envelope for dredge-up to take place. The previ- 
ous calibration of these quantities (Marigo et al. 1999) 
on the observed luminosity functions of C-stars has 
yielded (A = 0.50, logT$ led = 6.4) for the LMC and 
(A = 0.65, logT b drcd = 6.4) for the SMC. 

— Hot-bottom burning occurring in the most massive 
AGB stars (M > 4.5 M Q ), is followed in detail with 
the aid of complete envelope integrations (Marigo et 
al. 1998, Marigo 1998). The possible break-down of 
the core mass-luminosity relation and the CNO-cycle 
nucleosynthesis (possibly preventing the formation of 
luminous C-stars) are both taken into account. 

— Mass loss is included according to Vassiliadis & Wood's 
(1993) formalism, based on observations of pulsating 
AGB stars. In our calculations we assume that pulsa- 
tion occurs with periods (P) corresponding to either 
the fundamental mode (Po), or the first overtone mode 
(-Pi)- 

— Of great importance for the analysis developed in this 
work are the opacity prescriptions to be used in the 
envelope model. At high temperatures (T > 10000 K) 
we adopt Iglesias & Rogers (1993) for O- and C-rich 
mixtures. At low temperatures (T < 10 000 K) we use 
either the solar-scaled opacity tables by Alexander & 
Ferguson (1994; hereafter also Kfi x ), or the routine de- 
veloped by Marigo (2002; hereafter also K va r) that pre- 
dicts molecular opacities for any chemical composition 
of the gas. It should be remarked that before Marigo 
(2002), all TP-AGB tracks by Marigo et al. were cal- 
culated with Kfi x . 

In summary, for the purpose of our present work we 
will make use of two groups of synthetic TP-AGB calcu- 
lations, namely: 

— Old TP-AGB models with K fix (Marigo et al. 1999); 

— New TP-AGB models with K var (Marigo 2002). 

For the remaining prescriptions, we refer to the aforemen- 
tioned model outline unless otherwise specified. 

2.2. The formation of C-stars in other models and the 
role of molecular opacities 

Before proceeding with our analysis it is worth recalling 
briefly the current status of AGB models available in the 
literature and the importance of a consistent choice in the 
use of molecular opacities. 

First of all, the lack of efficient dredge-up in low-mass 
AGB models (say with M < 2M Q ) has been a diffi- 
culty shared by all full stellar evolution codes for a long 
time (the so-called "carbon star mystery" , as referred to 
by Iben (1981)). More recently this difficulty seems to 
be overcome by a few groups (e.g. Herwig et al. 1997; 
Straniero et al. 1997; Karakas et al. 2002), due to dif- 
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Fig. 3. Distribution of point sources from the 2MASS 
Second Incremental Data Release in the LMC region of 
the sky, projected on the (a, S) plane. To reduce confu- 
sion, we plot only the objects with K s < 12.5. The circle 
delimits the area selected for our analysis. 

ferent numerical treatments in describing the convective 
boundaries. 

However, even if some progress is achieved by full stel- 
lar evolution models in relation to the formation of low- 
mass C-stars, the distance between theory and observa- 
tions is still large. Currently there is no complete set of 
evolutionary tracks - based on full stellar models - ex- 
tending up to the end of the AGB phase for a sufficient 
coverage of stellar masses and metallicities. This prevents 
a systematic test of stellar models by comparison with 
observations, e.g. the luminosity functions of M- and C- 
stars, and simulations of CMDs including the most evolved 
AGB stars. The main reason is the heavy computational 
effort required to fully model the AGB phase. On the other 
hand, this limit can be totally removed if one allows for a 
certain loss of detail and opts for the more agile synthetic 
approach (see e.g. Groenewegen & de Jong 1993, Marigo 
et al. 1996, 1998 for recent works). 

Another heavy inadequacy of present stellar models 
(both full and synthetic) that has been recently high- 
lighted by Marigo (2002) is the use of fixed solar-scaled 
molecular opacities in AGB models experiencing the third 
dredge-up. 

In fact, this prescription (previously defined as ks x 
case) is still the standard choice common to most full 
and synthetic evolutionary calculations 1 of the AGB phase 

1 Actually, non-solar scaled molecular opacities were already 
taken into account in a simple way by Vassiliadis & Wood 
(1993), but no significant difference emerged in comparison 
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Fig. 4. 2MASS data for the innermost 4° of the LMC. 
Note the red tail of visible C-stars extending up to (J — 
K s ) ~ 2. 

(e.g. Mouhcine & Lancon 2002; Chieffi et al. 2001; Herwig 
2000; Ventura et al. 1999; Wagenhuber & Groenewegen 
1998; Forestini & Charbonnel 1997; Straniero et al. 1997). 

Marigo (2002) has proved, instead, that properly cou- 
pling the molecular opacities to the actual surface chem- 
ical abundances (K V ar) brings along such important con- 
sequences that the standard evolutionary scenario for the 
AGB evolution may be significantly affected. The reader 
is referred to that work for a detailed analysis. 

One major effect is that the surface enrichment in car- 
bon that leads to C/0> 1 - i.e. the formation of carbon 
stars - causes a significant increase in opacities, which in 
turn is responsible for a sudden and marked cooling of the 
stellar tracks in the H-R diagram. The sizeable excursion 
towards lower T e gt, displayed by the K var C-star tracks, 
sets the first significant difference compared to Kfi x C-star 
models that, instead, do not show any change in the slope 
of their evolutionary tracks while climbing the AGB at 
increasing luminosities. 

A closer comparison between models with Kfj x and 
K var , and the impact of the new opacities prescription will 
be discussed in Sect. 14.21 

3. The K vs. (J — K) diagram towards the LMC 

In this section, we aim to model the red tail of carbon 
stars in the LMC by using the TP-AGB models at our 
disposal. Our study is limited to a particular CMD - the 
K vs. (J — K) one - that reveals the red tail in all its 

to models with solar-scaled opacities, as the authors did not 
practically find carbon dredge-up in their AGB calculations. 
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prominence (see Fig.^). Needless to say, the same feature 
is displayed in a variety of different CMDs involving near- 
infrared passbands, as can be seen e.g. in van der Marel 
& Cioni (2001, their figure 2). 

3.1. The selected data 

In the following, we base our comparisons on 2MASS data 
towards the LMC, that present an extremely good photo- 
metric quality. More specifically, we use the point-source 
catalog of the 2MASS Second Incremental Data Release 
(see Cutri et al. 2002), which includes 10cr detections as 
faint as 16.3, 15.3, 14.7 mag in J, H and K s , respectively. 
As it can be seen in Fig. 03 these data wholly cover the 
LMC except for two limited strips of right ascension. In or- 
der to reduce the fore-/back-ground of Galactic stars and 
distant galaxies, we select only the stars located closer 
than 4 degrees (projected) to the LMC bar optical cen- 
ter (at a = 5 h 24 m , S = -69°44', cf. Nikolaev & Weinberg 
2000). This region, also considering the empty strip, trans- 
lates into to a total sky area of about 48.2 squared degrees. 
The corresponding CMD is shown in Fig. 

3.2. Simulations 

The K vs. (J — K) diagram has been simulated by using 
a population synthesis code (Girardi et al., in prepara- 
tion). In short, the code randomly generates stars follow- 
ing a given star formation rate (SFR), age-metallicity rela- 
tion (AMR) and initial mass function (IMF). The stellar 
intrinsic properties (luminosity L, effective temperature 
T e ff, surface gravity g, etc.) are interpolated over a large 
grid of stellar evolutionary tracks, based on Bertelli et al. 
(1994) for massive stars, Girardi et al. (2000) for low- and 
intermediate-mass stars, and complemented with the TP- 
AGB tracks from either Marigo et al. (1999) with Kfi x , 
or new calculations with K var made on purpose (refer to 
Marigo 2002; see also Sect.EJ. 

The evolution along the TP-AGB phase is initially de- 
scribed in terms of the properties in the quiescent stages 
between thermal pulses. Actually, there are two additional 
factors that cause significant excursions in luminosity and 
effective temperature from their quiescent values, namely: 

— occurrence of thermal pulses; 

— instability against pulsation. 

The former effect is included in our simulations. In 
practice, the luminosity variations driven by a thermal 
pulse (both the fast luminosity peak and long-lasting low- 
luminosity dip) are described according to the luminosity 
distributions over a pulse cycle as derived from Boothroyd 
& Sackmann (1988), and assuming that the envelope mass 
is the main factor determining the shape of the luminosity 
dip. Once a star is randomly scattered in logL according 
to this distribution, the corresponding excursion in effec- 
tive temperature is determined from a grid of envelope 
integrations performed in precedence. In this way, the in- 
trinsic dispersion of TP-AGB stars in the HR diagram is 



realistically simulated with just a modest computational 
cost. We recall that this dispersion is important mainly 
for low-mass AGB stars, with a maximum displacement 
of w —0.4 in log!/ (or ~ +1 mag), affecting up to 30 
percent of their TP-AGB lifetimes. 

As for the latter point related to stellar pulsation, it is 
worth premising the following. We recall that both 2MASS 
and DENIS provide single-epoch observations. Hence, the 
sampling of variable stars (e.g. Miras, SR variables) at a 
single and random phase of their variability cycle causes 
a further dispersion of data points in the observed CMD. 
This effect has not been simulated here for two reasons. 
First, we do not know exactly how L and T g vary during 
the pulsation cycle. Second, in any case the effect of this 
cyclical variation is easy to foresee: the bulk of TP-AGB 
stars will be scattered in the CMDs, to both higher and 
lower magnitudes and colours. The expected result is a 
sort of general blurring of the TP-AGB sequences in the 
CMD. Mira and SR variables have pulsation amplitudes 
ranging from 0.4 to 1.0 mag in the /-band (Hughes & 
Wood 1990). Therefore, to take this effect into account 
the simulated TP-AGB sequences should be blurred by 
just some tenths of magnitude. 

Once the stellar intrinsic properties are singled out, the 
photometry is generated by applying the extended tables 
of bolometric corrections from Girardi et al. (2002) for 
O-rich stars (with C/0< 1), and empirical relations - to 
be discussed in Sect. 14. ll below - for C-type stars (with 
C/0>1). 

Our simulations are computed in the J UK filter sys- 
tem defined by Bessell & Brett (1988). The expected dif- 
ferences from the 2MASS and DENIS systems - that use 
a "if-short" filter, K s - are probably very low, of less than 
say 0.1 mag in (J — K). This is illustrated in the case of 
2MASS by Cutri et al. (2002, section II,2d). 

We simulate the photometric errors by means of a 
Monte Carlo approach, adopting the distribution of er- 
ror values tabulated from 2MASS. We recall that tr errors 
are lower than 0.03 mag for K s Ss 13. 

Finally, the population synthesis code generates sim- 
ulations for both the LMC galaxy and the Galactic fore- 
ground, by using different distributions of reddening and 
distances. Results are shown in Fig. El for the specific case 
of TP-AGB models computed with solar-scaled molecular 
opacities (regx), and discussed in the following. 

3.2.1. The Galactic foreground 

The field foreground stars have been included with the 
aid of a complete Galaxy model (Girardi et al., in prepa- 
ration; Groenewegen et al. 2002). It consists of disk and 
halo components: 

The disc is described by a double-exponential density 
law, with a radial scale length (in Galactocentric distance) 
of 2.8 kpc, and a scale height H that increases with stellar 
age t as 

H = z Q (l+t/t ) a . 
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Fig. 5. Simulated CMDs for both the Galaxy foreground and the LMC population. Panel a: The Galaxy foreground 
stars. The coloured dots mark the stars belonging to the disk (green) and halo (black) populations. Panel b: The 
LMC population, marked in blue for all stars before the TP-AGB phase, and separated into O-rich (cyan) and C-rich 
stars (red) during the TP-AGB phase. Note that the C-star population distribute almost vertically over the red giant 
branch. 




J K 



Fig. 6. The same data as in Fig. EK, but in the Mr- vs. 
(J — K) plane. 



We adopt the parameters z = 95 pc, to = 0.5 Gyr and 
a = 1.66, that well describe the disk components as de- 
rived by Ng et al. (1997). It is worth noticing that this 
prescription allows for a significant increase of H at large 
ages, which is similar to assuming the presence of an old 



thick disk component. The distributions in mass, age and 
metallicity of the simulated disk stars are derived under 
the assumptions of: a constant star formation rate between 

12 Gyr ago and now; the empirical age- metallicity relation 
from Rocha-Pinto et al. (2000) with the addition of a dis- 
persion of 0.2 dex in [Fe/H] at any age; and the log-normal 
IMF from Chabricr (2001). 

The halo is described by an oblate spheroidal (cf. 
Gilmore 1984) of axial ratio q = 0.65, and a core radius 
of 2.8 Kpc. Halo stars are assumed to be old (from 12 to 

13 Gyr), and to follow a Gaussian distribution of metal- 
licities of mean [Fe/H] = — 1.6 and standard deviation of 
1.0 dex. 

This set of prescriptions has been constrained by 
Groenewegen et al. (2002), by using a set of Deep 
Multicolor Survey (Osmer et al. 1998) and ESO Imaging 
Survey data (Prandoni et al. 1999). Interestingly, most of 
the calibration data were in visible pass-bands, whereas 
the model is now being applied to the interpretation of 
near-infrared data. 

For the purpose of our study, we simulate a conic sec- 
tion of the Galaxy occupying the same sky area as in the 
selected 2MASS data. The field center has galactic coor- 
dinates £ = 280°.46, b = -32°.89. We recall that, for this 
particular line-of-sight, other Galactic components (spiral 
arms, the bulge) are of no relevance. 

An important aspect is that in the modelling we can 
distinguish different kinds of stars and hence tackle the ori- 
gin of the several observed features in the CMD of 2MASS. 
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Our field simulation is presented in (Fig. . In this 
plot, we can identify three marked vertical sequences of 
foreground disk stars (green dots), located at (J — K) 
colours of about 0.35, 0.65, and 0.9. Their origin is easily 
understood by looking at the corresponding plot in ab- 
solute magnitude (Fig. |BJ: The sequence to the left (at 
(J — K) ~ 0.35) is defined mainly by the old disk turn-off, 
with typical masses of M ~ 0.9 M ; the central one (at 
(J — K) ~ 0.65) hosts fainter RGB and red clump stars; 
and the right one (at (J — K) ~ 0.9) corresponds to low- 
mass dwarfs with M Ss 0.6 M Q . The vertical development 
of these strips in the CMD (Fig. [SJi) simply reflects the 
large range of distances involved. 

We notice that disk stars younger than ~ 4 Gyr (M 
1.3 Mq) do not populate significantly this diagram, since 
they are found preferentially at small scale heights and 
hence at lower galactic latitudes and brighter apparent 
magnitudes (i.e. K <5 10). The contribution of halo stars 
(black points) is also of minor importance, accounting for 
just a handful of points in the synthetic CMD, mostly at 
fainter apparent magnitudes (i.e. K ^ 13), due to their 
large mean distances. 

Finally, an important point is that field stars are not 
expected to reach colours redder than (J — K) ~ 1.0, so 
that they hardly contaminate the CMD features produced 
by the LMC population of AGB stars. Their intepretation 
is our next goal. 

3.2.2. The LMC population 

The LMC population has been initially described by 
means of a very simplified model, under the assumptions 
of i) a fixed distance of 52.3 Kpc (Girardi & Salaris 2001), 
ii) constant SFR from 15 Gyr ago to now, and iii) a con- 
stant metallicity, Z = 0.008, at any age. The latter choice 
is made for the sake of simplicity, since it requires the use 
of a single set of evolutionary tracks, hence relieving us, 
for the moment, from the effort to compute complete sets 
of TP- AGB models for several metallicities 2 . Anyway, the 
quality of our results should not depend much on these 
assumptions 3 . 

In order to get a number of LMC giants comparable 
to the observed sample (Fig. 0J, we have to simulate a 
total star formation of 5 x 10 8 M Q during the entire LMC 
history. Of course, this is just a rough estimate of the 
mass consumed to form stars, its value depending on the 
adopted assumptions about SFR (constant), IMF (from 
Kroupa 2001), and size of the simulated region (i.e. within 
the innermost 4° of the LMC field). Anyway, it is inter- 

2 Extended sets of TP- AGB tracks and isochrones for a wide 
range of metallicities are being computed and are soon to be- 
come available. 

3 In one of our test cases, we used TP- AGB tracks interpo- 
lated within an extended grid of 5 different metallicities, to- 
gether with a more realistic AMR (from Pagel & Tautvaisiene 
1998). No significant difference was noticed in the resulting 
CMDs. 



esting to notice that this number is just an order of mag- 
nitude lower than the total LMC mass (including dark 
matter; see van der Marel et al. 2002), and comparable to 
the present LMC Hi mass (Staveley-Smith et al. 2003). 

The simulated LMC population draws a few almost 
vertical features in the CMD (Fig. EJs), mainly popu- 
lated by: young and intermediate-age main-sequence stars 
(with (J — K) < 0.5), intermediate-mass stars in the 
stage of core-He burning (the long diagonal sequence at 
(J — A") ~ 0.8), intermediate-age and old RGB stars of 
larger luminosities (the prominent sequence at (J—K) ~ 1 
and with K > 12), and intermediate-mass stars during the 
early- AGB phase (the weaker sequence departing slightly 
from the left of the RGB-tip, with the same inclination 
as the RGB). All these features (marked with blue dots) 
correspond to O-rich stars in evolutionary stages previous 
to the TP- AGB. 

The O-rich TP-AGB stars (marked with cyan dots) 
trace two well-defined sequences: the first consists of 
intermediate-mass stars that somewhat extend the plume 
of early- AGB stars towards higher luminosities; the second 
defines a clear strip located directly above the RGB-tip, 
corresponding to low-mass stars. Finally, The C-rich TP- 
AGB stars (marked with red dots) distribute above and 
to the right of these latter sequences. 

Combining both panels of Fig. it turns out that 
our initial simulations explain quite well the several "fin- 
gers" displayed by the 2MASS data of Fig. fj] but for one 
major discrepancy: the location of C-stars, predicted by 
these particular models, simply align above the sequence 
of LMC O-rich giants, failing to form the red tail. As for 
the other CMD features, the agreement with data is indeed 
quite good, and our interpretations substantially agree 
with those advanced by Nikolaev & Weinberg (2000). 

4. Probing the red tail of C-stars 

In the following we will focus on the region of our sim- 
ulated diagrams populated by TP-AGB stars, with the 
aim of casting light on the main cause that produces the 
observed red tail of C-stars in near-infrared colors. We 
investigate two possible factors separately, namely: the 
colour-T e ff relations (derived from empirical calibrations), 
and the evolution in effective temperatures (derived from 
stellar models). 

4.1. Testing different colour-T e s relations 

Let us start by analysing the effect on our simulations of 
different (J — K) — T e g transformations. It is important to 
remark that all computations discussed in this section are 
carried out by adopting TP-AGB tracks for solar-scaled 
Kfi x opacities. Results are shown in Fig. [7| 

As a first step, in Sect. 14.1.11 we consider the sharp 
dichotomy between colour-T^g- relations derived from ob- 
served M-type stars (e.g. Fluks et al. 1994) and C-stars 
(e.g. Bergeat et al. 2001). Such differences essentially re- 
flect the abrupt change in the dominant molecular species 
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- producing the line blanketing - when passing from an 
O-rich spectrum (e.g. TiO, VO, H2O) to a C-rich one (e.g. 
C a , CN, SiC) 

Then, within the class of C-stars, in Sect. 14.1.21 we 
attempt to account, in a simple way, for the possible com- 
position dependence (related to the extent of carbon en- 
richment) of colour-Teff relations. 

4.1.1. A T eff - (J - K) relation valid for Galactic 
C-stars 

Figure^ shows our initial simulation: The TP-AGB mod- 
els are calculated with K-fix molecular opacities, and the 
conversion from T e ff to (J — K) is based on Fluks et al. 
(1994) empirical spectra and temperature scale for nearby 
M giants (see Girardi et al. 2002 for details). This means 
that a relation valid for O-rich stars has been applied also 
to C-stars. The discrepancy with observed data is obvious: 
C-stars just align above the sequence of O-rich stars, and 
do not develop any hint of red tail. 

As a first attempt to include proper relations valid 
for C-stars, we apply the Bergeat et al. (2001) formulas, 
that give -ftf-band bolometric corrections and (J — K) as 
a function of T e g. These fitting relations are derived from 
a sample of Galactic C-stars. In this context we would 
like to note that the Bergeat et al. (2001) relationship is 
essentially the same as that of Bessell et al. (1981). The 
resulting simulations are shown in Fig. |7p. With respect 
to the previous case, the sequence of C-stars is shifted to 
the blue, and not to the red as required to explain the 
data. 

In this and other plots we show next, we have also 
adopted the empirical X-band bolometric corrections for 
C stars presented by Frogel et al. (1990). They differ lit- 
tle (by ~ 0.1 mag) from those independently derived by 
Bergeat et al. (2001), and yield closely similar simulated 
CMDs. 

4.1.2. The effect of molecular blanketing for C-stars 

Molecular blanketing is responsible for the distinction be- 
tween M- and C-spectral features and also affects the in- 
frared colours of C-stars themselves (e.g. Cohen et al. 
1981; Loidl et al. 2001). In this latter case, the extent of 
molecular blanketing is essentially related to the amount 
of carbon available to form molecules other than CO. 

In their study of Galactic and LMC C-stars, Cohen et 
al. (1981) pointed out that the observed mean (J — K) 
colour increases steadily as the carbon abundance class - 
defined by the C2 band strength - increases (from say C,2 
to C,5). Moreover, the authors showed that, for a given ef- 
fective temperature, redder JHK colours are predicted by 
model atmospheres with larger C2 and CN enhancements. 

In our work we attempt to account, in a simple way, 
for the effect of molecules on the infrared colours of C- 
stars, by introducing a C/O-dependence into the empirical 
(J — K) — T off relation by Bergeat et al. (2001). The C/O 



data for Galactic C-stars are taken from Lambert et al. 
(1986). Figure [3] displays the observed trends of the (J — 
K) colour as a function of T e ff and C/O ratio. The linear 
bi-parametric fit relation to the data is: 

(J — K) = 17.32 - 4.56 logT cff + 0.052 C/O (1) 

We see that the (J — K) — T e g relation (top-left panel of 
Fig. |SJ is quite narrow for Galactic C-stars, and it is well 
reproduced by our fitting formula for C/0= 1.1, in agree- 
ment with the observational finding that most Galactic 
carbon stars have low C/O ratios, not much larger than 
one (Lambert ct al. 1986). Moreover, despite some ob- 
served scatter and the small number of data points, we 
also note that there is a positive correlation between the 
(J — K) colour and the C/O ratio (bottom-left panel of 
Fig. |SJ . This trend is included in our fitting relation as it 
predicts a systematic shift towards redder (J — K) colours 
for e.g. C/0= 3.0, compared to the case C/0= 1.1 (top 
left panel). 

Overall, Eq. fits the data to better than 0.1 mag in 
(J — K), as indicated by the residuals in the right panels 
of Fig.il 

Actually, in a recent study Matsuura et al. (2002) 
have pointed out that the molecular features in the spec- 
tra of LMC C-stars are consistent with C/O ratios larger 
than those of Galactic C-stars. This empirical finding com- 
pletely supports the results of our TP-AGB calculations, 
as shown in Fig. EI The same figure indicates that, in or- 
der to model the LMC population, we will have to apply 
Eq. (0) to C/O values as high as 4, that is beyond the 
maximum observed values of Galactic C-stars. 

Thus, by using Eq. Q in our simulations, we get the 
synthetic CMD of Fig. 0;. As it can be noticed, the de- 
pendence of (J — K) on the C/O ratio now causes the 
sequence of C-stars to be significantly redder, reaching 
(J — K) = 1.7. However, the C-stars are still largely dis- 
tributed among and above the O-rich giants, and no clear 
red tail is drawn. 

If we limit consideration to models with C/0< 2, i.e. 
plotting only those stars for which Eq. is strictly valid, 
we get Fig. [7|i. In this case, we can notice that the low- 
C/O C-stars contaminate significantly the CMD region 
corresponding to O-rich AGB stars. On the contrary, such 
a high fraction of C-stars mixed up with M-stars is not 
present in observed samples. It is clear that this discrepant 
feature would remain for whatever hypothetical (J — K) — 
T e ff— C/O that could push somehow the models to reach 
(J — K) colours as high as 2.0. 

In conclusion, our test calculations seem to indicate 
that the principal cause of the failure of models with solar- 
scaled molecular opacities in reproducing the red tail of 
carbon stars should not be ascribed to the (J — K) — T e ff 
relation. 

4.2. Testing different TP-AGB tracks 

If we assume that Eq. reasonably well describes the 
(J — K) colours of C-stars, as supported by the low fitting 
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Fig. 7. Simulated CMDs based on TP-AGB tracks with k^ x prescription for molecular opacities. The simulations 
includes both the Galaxy foreground and the LMC population. The coloured dots mark several kinds of stars: fore- 
ground stars in the Galaxy disk (green) and halo (black) ; LMC giants before the TP-AGB phase (blue) , and separated 
into O-rich (cyan) and C-rich stars (red) during the TP-AGB phase. Note that the C-star population distribute almost 
vertically over the red giant branch. Different (J — K) vs. T c g transformations are tested in the simulations of the C 
stars: Panel a: The same as for M stars (based on Fluks et al. 1994); Panel b: Bergeat et al. (2001) one; Panel c: 
A fitting relation involving C/O (Eq.QJ. A modest inclination appears in the C-star sequence. Panel d: The same as 
in Panel C, but plotting only the data with a maximum C/O value of 2.0. 



residuals (of ^0.1 mag, see the right panels of Fig.[BJl, we 
are left with just one alternative: redder (J — K) colours 
could be attained only by assigning lower T c g to C-star 
models. Moreover, in order to locate the C-star models 
onto a different branch from that of the O-rich models, a 
net T e g segregation between the two classes turns out to 



be necessary. These requirements are indeed fulfilled by 
Marigo (2002) TP-AGB models with variable molecular 
opacities. 

Figure ^3 illustrates a set of Z = 0.008 isochrones 
based on Marigo's (2002) models. Similar isochrone sets 
have been computed for a large variety of dredge-up pa- 
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Fig. 8. Left panels: (J — X) vs. T e g and (J — K) vs. C/O relations for Galactic C-stars. Observed data (filled 
squares) are from Bergeat et al. (2001, objects common to their tables 4 and 13). In the upper panel, lines show the 
predictions of the fit relation given by Eq. QJ, assuming a constant C/O ratio, equal to 1.1 (solid) and 3.0 (dashed). 
Right panels: The residuals of fitting Eq. JIJ to the data, plotted as a function of both T c g and C/O. The data is 
well reproduced by this relation to within 0.1 magnitudes. 



rameters and metallicities, and will be extensively de- 
scribed in a future paper. The conversion to (J — K) 
colours for C-star models is now based on Eq. (JJJ . It is also 
interesting to notice that in this case we could not even 
have applied the same colour transformations derived for 
M-stars, since the effective temperatures of C-star mod- 
els with already much cooler than the validity 
lower- limits of such relations. 

4.2.1. The sole effect of variable molecular opacities 

The effect of TP-AGB models with variable molecular 
opacities is illustrated in the simulations of Fig. ^2 TP- 
AGB models are computed with different assumptions (see 
also Sect. regarding the onset of the third dredge-up - 
controlled by the parameter T dred -, and the onset of the 
super-wind regime - depending on the pulsation mode, ei- 
ther fundamental (P = Pq), or first overtone (P = Pi). In 
all cases, the sequence of C-stars clearly departs from the 
branch of O-rich stars. 

In these K var models the excursion towards redder 
(J — K) colours is mostly caused by the displacement 
towards cooler T e g. Instead, the increase of the C/O ra- 
tio during the TP-AGB evolution should play a smaller 
role. In fact, considering that the maximum C/O ratio 
reached by these models is of about 4 (Fig. (SJ, the C/O- 
dependence in the colour transformation of Eq. ac- 
counts for at most 0.2 mag of the total colour excursion 
at the end of the C-star evolution, the mean value being 
of about 0.1 mag. 



In Fig. ITTk . we adopt the same dredge-up parameters 
as calibrated in Marigo et al. (1999), i.e. logT^ rcd = 6.4 
and A = 0.5. The simulated C-star red tail resembles very 
much the observed one, but for two aspects. 

First, the red tail shows up at a too low luminosity. 
It starts at K = 11.6 when (J — K) = 1.2, whereas the 
observed tail is at least 0.6 mag brighter. This problem 
can be easily solved with a re-calibration of the dredge-up 
parameters. To give an example, Fig. lllb shows a simula- 
tion that makes use of TP-AGB models recomputed with 
logT 1 drcd = 6.45. This parameter choice has the effect to 
delay the onset of the third dredge-up to higher stellar 
core masses and hence higher luminosities. As we can see, 
now the red tail starts at K = 11.0 when (J — K) = 1.2, 
which is about the right luminosity level to fit the data. 

It is worth remarking that in both simulations pre- 
sented in the upper panels (a and b) of Fig. 1111 we have 
assumed that all TP-AGB stars are fundamental-mode 
pulsators, meaning that the Vassiliadis & Wood's (1993) 
mass-loss formalism is applied with M(P) = M(Pq). This 
assumption introduces the second problem of the models, 
which is discussed below. 

4.2.2. Extending the red tail with first-overtone 
pulsators 

The mentioned second theoretical difficulty resides in 
the extension of the red tail, which in both simulations 
(Figs. II lh and b) clearly reaches (J — K) — 1.8, but 
hardly (J — K) = 2.0 as in the observations. Although 
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Fig. 11. Simulated H-R diagram based on TP-AGB tracks with variable molecular opacities (i.e. Marigo 2002). Note 
the departure of the C-star population away from the red giant branch towards redder colours. The mass loss pre- 
scriptions is taken from Vassiliadis & Wood (1993), assuming that AGB stars are either fundamental-mode or first 
overtone pulsators. Panel a: case logT drcd = 6.4, fundamental mode. Panel b: case logT' drcd = 6.45, fundamental 
mode. Panel c: case logT drcd = 6.4, first overtone. Panel d: case logT dred = 6.45, first overtone. 



this point could be considered as a minor discrepancy (at 
least compared with the Kfi x cases illustrated in Fig. 0) , it 
deserves some attention. There are at least two possible 
ways out. The first is that we have not considered the effect 
of circumstellar reddening on the colours of C-stars. This 
should be ascribed to thick circumstellar material ejected 
by the coolest C-stars during an extreme mass-loss phase. 
In fact, a few LMC C-stars are known to extend up to 
(J - K) - 6 (Cioni ct al. 1999; Nikolaev & Weinberg 



2000); they are among the so-called "obscured" or "dust- 
enshrouded" C-stars identified with IRAS sources (see van 
Loon 1999), which are found to undergo heavy mass- loss. 
The question, however, is if this sort of self- reddening is al- 
ready effective in red tail C-stars with say (J — K ) ^ 1.8, 
being able to shift their colours up to (J — K) = 2.0. This 
point has to be investigated with the aid of more detailed 
modelling. 
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Fig. 9. Final C/0 ratio (at the end of the TP-AGB phase) 
as a function of the stellar mass (at the onset of the TP- 
AGB phase), predicted by synthetic models with variable 
molecular opacities K var , and parameter set (A = 0.50, 
logT drod = g 4 ) Results for solar metallicity (Z = 0.019) 

are compared to those for LMC composition (Z = 0.008). 
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Fig. 10. Z = 0.008 theoretical isochrones for TP-AGB 
models computed with variable molecular opacities. The 
dredge-up parameters are the same as in the models of 
Fig. The dotted lines mark the isochrone sections cor- 
responding to C-stars. The conversion to the (J — K) pho- 
tometry is now based on Eq. 



The second explanation stands on the assumption that 
the dominant pulsation mode of AGB stars is the first 
overtone (with period P = Pi). This would allow us to 
get simulated colours further increased by 0.2 mag, while 
keeping the same T e g vs. (J — K) relation. This is illus- 
trated in the bottom panels of Fig. ^] For these simu- 
lations we use TP-AGB models computed with the same 
dredge-up parameters as in the upper panels, but assum- 
ing that Mira and SR variables are first overtone pulsators. 
Under this hypothesis, variable AGB stars should have 
shorter pulsation periods for given stellar L and T e g, com- 
pared to the fundamental- mode counterparts. 

As a consequence, by virtue of the empirical posi- 
tive correlation between the mass-loss rates and the pe- 
riods of Mira stars, we expect that for the shorter first 
overtone periods (Pq/Pi ~ 2.2) the onset of the super- 
wind regime is postponed, thus implying that the end of 
the AGB evolution is delayed. We quantify this effect by 
adopting the Vassiliadis & Wood's (1993) formalism with 
M(P) = M(Pi) (Figs. HTb.dV In this way we find that 
C-star models are allowed to reach slightly higher L and 
lower T c ff, and their lifetimes increase accordingly. This 
latter effect is also relevant if we consider that longer life- 
times are needed to better reproduce the observed C-star 
counts in LMC clusters (Girardi & Marigo 2003, and work 
in preparation). 

In Fig.lllb. we can notice that now the red tail actually 
extends up to (J — K) = 2.0. However, it still appears at 
too faint luminosities, due to the low value assumed for 

Instead, for a larger 



7^ red (see remarks in Sect. I4TT 
ydrcd^ consequent later onset of the third dredge-up is 
able to shift the entire red tail to the right luminosities, 
as shown in Fig. IllH . 

From a comparison with our reference plot of Fig. ^ 
we can conclude that the simulation of Fig. Illt l best re- 
produces the observed morphology of the red tail. On the 
other hand, the introduction of the first overtone hypoth- 
esis produces an undesirable effect in the sequence of old 
O-rich stars: they now extend to much higher luminosi- 
ties (up to K = 10.4 and K = 10.0 in Figs. HTJ: and d, 
respectively) than observed (up to K = 11.2 in Fig. 0}. 

As a matter of fact, nowadays the observational sce- 
nario of long-period variables (LPVs) appears quite com- 
plex, mainly due to the observed mixture of different pul- 
sation modes, and the lively controversy about their as- 
signment (see e.g. Wood et al. 1999, Feast 1999, Whitelock 
& Feast 2000, Cioni et al. 2001, Noda et al. 2002). The 
situation gets even more intricate if one considers the ev- 
idence that some LPVs are found to be multi-period pul- 
sators (Bergeat et al. 2002), and the possibility of an evo- 
lutionary path across different period-luminosity relations, 
so that LPVs may switch pulsation modes while evolving 
along the AGB (Cioni et al. 2001). 

Without entering the open debate on the dominant 
pulsation mode, we perform a very simple test on the basis 
of purely theoretical considerations. Considering that the 
red tail of C-stars is somewhat better reproduced with the 
first-overtone assumption which, instead, does not seem to 
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suitably describe the O-rich branch, we assume a mixed 
population of fundamental-mode and first-overtone pul- 
sators, with the following composition: 

1. All stars with mass lower than 1.3 Mq, and larger than 
3.5 Mq, are fundamental-mode pulsators just because 
they do not evolve into a C-star phase; 

2. The remaining stars are assigned either as 
fundamental- mode or as first-overtone pulsators. 
We assume a mix of 50 percent between the two types. 

Assumption J5J may be not realistic enough, because the 
simulated stars do not switch between pulsation modes 
as they evolve, but they can just follow one between two 
different evolutionary paths. If this is an oversimplifica- 
tion of the problem - unavoidable because the mechanisms 
leading to change the pulsation mode cannot be presently 
modelled - at the same time it is a resonable and instruc- 
tive experiment. In Fig. I12l we present the resulting simu- 
lation, compared with the selected 2MASS data. We can 
notice, in fact, that the resulting sequence of O-rich giants 
is noticeably shortened in luminosity, as expected. 

We may therefore conclude that now our final simu- 
lation actually meets 2MASS and DENIS observations, 
as all CMD features are well reproduced by the models, 
including the C-star red tail and the plume of O-rich low- 
mass TP-AGB stars. 

5. Concluding remarks 

In this paper, we have shown that most features seen in the 
2MASS and DENIS K s vs. (J-K s ) diagram are quite well 
described by our present stellar models. In particular, the 
observed tail towards the reddest (J—K s ) colours provides 
a direct and stringent probe to the evolutionary theory of 
TP-AGB stars. In this regard, we have pointed out that: 

1. The red tail of carbon stars is missed by models in 
which the TP-AGB evolution is followed using solar- 
scaled molecular opacities. 

2. This failure can hardly be caused by some hypothetical 
inadequacy in the present-day T e g vs. (J—K) relations 
for C-stars. Our tests indicate that in order to reach 
the observed (J — K) colours of C-stars with "solar- 
scaled opacity" models, quite high mean values of the 
C/O ratio should be exhibited by the bulk of carbon 
stars, which seems unrealistical. 

3. Instead, the red tail comes out naturally of TP-AGB 
models which are calculated with variable molecular 
opacities (cf. Marigo 2002), i.e. opacities consistently 
computed for the current chemical composition of the 
stellar envelope. This can be understood considering 
that, as the third dredge-up makes the C/O ratio in- 
crease above 1, drastic changes in the main sources of 
molecular opacities cause both the transition from M- 
to C-spectral types, and a significant decrease of the 
effective temperature. The latter aspect has been so 
far largely ignored in AGB evolution models. Contrary 



to what claimed by Mouhcine & Langon (2002), the 
present study demonstrates that in order to model 
properly the properties of AGB, in particular C-rich, 
stellar populations, one cannot ignore the large opac- 
ity effects have on the effective temperature of C-stars, 
and more generally on the AGB evolution. 

4. Our simulations with "variable opacity" models pro- 
vide a few additional useful indications about the TP- 
AGB evolution. For instance, models with logT^ rcd = 
6.45 seem to reproduce well the -RT-band luminosity 
of red tail C-stars. This means that the onset of the 
third dredge up should occur at slightly larger core 
masses, as compared to previous models by Marigo et 
al. (1999). 

5. Moreover, in the context of the present models, based 
on the use of Vassiliadis & Wood (1993) prescription 
for mass loss and an empirical T c g vs. (J — K) rela- 
tion, we may derive two indications about the pulsa- 
tion modes of variable AGB stars. First, the colour 
extension of the red tail may imply a significant frac- 
tion of first-overtone pulsators among C-type stars. 
However, the same extension could be reproduced with 
fundamental- mode pulsators if circumstellar reddening 
starts to be important for C-stars with (J — K)o ~ 1.8, 
a point that has still to be properly investigated. 
Second, the luminosity extension of the sequence popu- 
lated by O-rich low-mass TP-AGB stars would suggest 
that most of M-type variables are fundamental-mode 
pulsators. 

All these indications are being used to calibrate our 
TP-AGB models, from which extended sets of theoretical 
isochrones at varying metallicity will be derived and soon 
made available. Of course, new models will be computed 
for the variable-opacity case only. 
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